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Abstract 

In supersymmetric Grand Unified Theories, proton decay medi- 
ated by the color-triplet higgsino is generally problematic and requires 
some fine-tuning of parameters. We present a mechanism which nat- 
urally suppresses such dimension 5 operators in the context of SUSY 
SO (10). The mechanism, which implements natural doublet-triplet 
splitting using the adjoint higgs, converts these dimension 5 opera- 
tors effectively into dimension 6. By explicitly computing the higgs 
spectrum and the resulting threshold uncertainties we show that the 
successful prediction of sin 2 9\y is maintained as a prediction in this 
scheme. It is argued that only a weak suppression of the higgsino me- 
diated proton decay is achievable within SUSY SU (5) without fine- 
tuning, in contrast to a strong suppression in SUSY SO (10). 
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1. Introduction 



The dramatically precise unification of couplings that occurs in the minimal 
super symmetric SU(5) model has been much touted, and indeed is strik- 
ing. A fit 2 to all W,Z, and neutral current data (using m t = 138GeU and 
m H = M z ) gives sin 2 9 W (M Z ) = 0.2324 ± 0.0003, while in SUSY SU(5) one 
has 2 sin 2 Q W {M Z ) = 0.2334±0.0036 (we have combined the uncertainties due 
to a s (Mz), a(Mz), sparticle thresholds, m t ,mh , high-scale thresholds and 
non-renormalizable operators). As a measure of how significant this agree- 
ment is consider that the addition of just one extra pair of light higgs doublet 
supermultiplets would increase the SUSY GUT prediction of sin 2 Q W (M Z ) to 
about 0.256, leading to a discrepancy of over six standard deviations. 

At the same time SUSY SU (5) suffers from a problem 3 with proton decay 
arising from dimension-5 operators caused by the exchange of color-triplet 
higgsinos in the same SUSY-Sf/ (5) multiplet with the light higgs doublets, 
H and H'. For 'central values' of the model parameters the predicted proton 
life-time mediated by the dimension-5 operators 4 ' 5 is shorter than the current 
experimental limits. 6 A certain amount of adjustment is then required for 
consistency, which pushes parameters to the corner of their allowed region. 
For short we will henceforth refer to this as the 'higgsino-mediated proton 
decay (HMPD) problem'. 

It would seem that one cannot take seriously the unification of couplings 
of SUSY-GUTs, however impressive, in the absence of satisfactory mecha- 
nism that suppresses higgsino-mediated proton decay. Two requirements for 
a 'satisfactory' mechanism that we will impose are that it involves no 'fine- 
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tuning' or artificial adjustments of parameters, and that the unification of 
couplings is maintained as a prediction. (We emphasize that word because a 
discrepancy in sin 2 9w can often be remedied by introducing ad hoc new par- 
ticles, threshold effects, etc; but we would not regard the resulting agreement 
as being in any sense a prediction.) 

The higgsino-mediated proton decay problem 3 is easily described. In 
SUSY SU(5) models there exists a pair of higgs super-multiplets, that we 
will denote 5 H +5' H . Under the standard model group, Gs = SU (3) x SU (2) x 
U(l), these decompose as 5# = {(1,2,+|) + (3, 1, — |)} = {2 h + 3h}, and 
5' H = {(1, 2, -\) + (3, 1, +§)} = {2' H + 3^}. The 2 H and 2~' H are just the 
familiar H and H' of the supersymmetric standard model, and their cou- 
plings to the light quarks and leptons are therefore fixed in terms of their 
vacuum expectation values and the light fermion masses. Then by SU (5) the 
Yukawa couplings of the 3h and 3'^ are fixed as well. If there is a Dirac mass 
term connecting the higgsinos in 3# and 3' H to each other, i.e., a term of 
the form M(3h3' h ), then the diagram shown in figure 1 exists, which depicts 
a baryon-number violating process mediated by colored higgsino exchange. 
The higgsino, of course, converts the quarks and leptons to their scalar part- 
ners, so figure 1 needs to be 'dressed' for it to correspond to proton decay. 
Dressing by VU-ino exchange is by far the most dominant, the resulting life- 
time for p — > K + u fl for example (the anti-symmetry of the relevant operator 
causes the change in flavor) has been estimated 4 to be 
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T(p^K+v») = 6.9xl0 2 Vx|(^f^)(^)(^ff) 

( M Hc \ TeV- 1 2 
\l0 16 GeV) f(u,d)+f(u,e) ' 

Here f3 is the relevant nuclear matrix element which lies in the range f3 = 
(0.003 — 0.03)GeV 3 . A s is the short- distance renormalization factor (A s ~ 
0.6), tan/3# — (H) / (H'), and y tk parametrizes the contribution of the top 
family relative to the first two (0.1 <| y tk |< 1.3 for m t = WOGeV). The 
functions / arise from the one-loop integrals and f(u,d) ~ m^ ± /m| for 
m w± < rriQ. 

The prediction (1) is to be compared with the present experimental lower 
limit r(p — > K + u fJi ) > 1 x 10 32 yr. 6 There is already, as can be seen, somewhat 
of a difficulty reconciling these numbers. At least several of the following 
conditions should be fulfilled: (a) The nuclear matrix element (3 is near 
the lower end, (3 ~ 0.003 GeV 3 , (b) tan(3 H is not too large, (c) either the 
W ± is significantly lighter than the squark Q or vice versa, (d) the colored 
higgsino mass should exceed the GUT scale significantly (M# c > 10 l7 GeV), 
(e) there is some cancellation in (1 + y tk ) between the third family and the 
first two family contributions. Obviously this pushes almost all parameters 
to their corners. Although not excluded, the problem begs for a more elegant 
explanation. 

The necessity of the Dirac mass term M(3#3#) for obtaining the D = 5 
baryon-violating operator is a crucial point. A quartic term in the chiral 
superfields will be suppressed by M ^ UT or — depending on whether it is an 
F-term or a D-term. To get an F-term all the left-handed chiral superfields 



4 



must be coming into the diagram (or out of it). This requires the chirality- 
flipping mass insertion coming from the M(3h3' h ) term as shown in fig. 1. 
If the mass M vanished, one could only write superfield diagrams like fig. 
2, which clearly give D-terms and therefore are suppressed by -jA — . (The 

GUT 

F-terms correspond to higgsino exchange, the D-terms to higgs exchange.) 

The foregoing considerations have suggested to several authors 7 ' 8 an ap- 
proach to suppressing higgsino-mediated proton decay by imposing a (Peccei- 
Quinn type) symmetry that suppresses the Dirac mass between the 3h and 
3' H . However, this approach leads to a dilemma. In the simplest SUSY SU (5) 
model it is precisely the M(3h3' h ) Dirac term that gives to the 3h and 3'h 
the superlarge mass that they must have (else even the D = 6 operators 
would cause a disaster). How, then, can the 3# and 3' H be made superheavy 
and yet not have a large Dirac mass connecting them to each other? To 
resolve this, one ends up introducing new color-triplet higgs superfields, a 
3h to mate with the 3h and a 3'^ to mate with the 3' H , so that all color 
triplets become heavy while still leaving the unprimed and primed sectors 
disconnected. That means having four instead of the minimal two 5-plets of 
higgs. The situation can be represented diagrammatically as follows 




5 H 5h 5'h 5' h 
where the solid horizontal lines representing superlarge Dirac masses. The 2 
and 2' are the usual light H and H' of the MSSM. But now it is to be noticed 
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that there are two additional light doublets, 2 and 2'. As noted above, this is 
unacceptable if the dramatic unification of gauge couplings is to be preserved 
as a prediction. 

One could remove the extra pair of Higgs(ino) doublets to superlarge 
scales by introducing a mass term M(5h5' h ), indicated by the dashed line 
in eq. (2). This, of course, would reintroduce the higgsino-mediated proton 
decay as shown in fig. 3. The parameter M would then control both the 
higgsino-mediated proton decay and the mass of the pair of extra doublets. 

As we shall see in the later sections, this situation is typical of mech- 
anisms that naturally suppress higgsino-mediated proton decay: (a) there 
is a doubling of the higgs sector, (b) there is the consequent danger to the 
unification of couplings of extra light fields in incomplete multiplets, and (c) 
there is a parameter M which controls both proton decay and the mass of 
these extra fields. 

Two approaches, therefore, appear to be possible. (1) Weak suppression 
of higgsino-mediated proton decay, resulting from M being of order - but 
slightly less than - M GUT . For example, with M = jqM gut higgsino- 
mediated proton decay is suppressed by a factor of 1CT 2 while at the same 
time those extra fields whose mass is given by M will only lead to small 
threshold corrections to sin 2 9w- 111 this case the suppression of higgsino- 
mediated proton decay is just numerical; there is no symmetry or other 
qualitative explanation of it. One would have no a priori reason therefore to 
expect the suppression to be particularly large. A hope would therefore exist 
that p — > K + v^,n — > K°v, etc. might be seen at super-Kamiokande. (2) 
Strong suppression of higgsino-mediated proton decay would result if (due to 
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some approximate symmetry perhaps) M were much less than Mqut, say, 
0(M W ). In that case it is imperative that there be no 'extra' fields (i.e., 
beyond the minimal supersymmetric standard model) in incomplete SU(5) 
multiplets whose mass is proportional to M, or else the unification of cou- 
plings would be destroyed. To achieve this without fine-tuning turns out 
to be a non-trivial problem. One of the main conclusions of this paper is 
that such a strong suppression of proton decay can only be achieved in a 
satisfactory way in SO (10) (or larger groups). 

The whole problem of higgsino-mediated proton decay is of course inti- 
mately connected to the well known question of 'doublet-triplet' splitting. 9-13 
We will show that the most satisfactory treatments of this problem make use 
of an old but somewhat neglected idea for doublet-triplet splitting in SO (10) 
using the adjoint higgs due to Dimopoulos and Wilczek. 13 

Our paper is organized as follows. In section II we review the Dimopoulos- 
Wilczek idea for doublet-triplet splitting and show how both weak suppres- 
sion and strong suppression of HMPD can be achieved naturally in SUSY 
SO (10). In section III we consider SUSY SU(5) and show that only weak 
suppression of HMPD can be achieved without fine-tuning parameters. In 
section IV we discuss flipped SU(5). In section V a closer examination of 
the Dimopoulos- Wilczek mechanism in SO (10) is undertaken, and we show 
that it can be made viable and consistent. Our conclusions are summarized 
in section VI. In appendix A we give the details of the minimization of a 
realistic SO (10) superpotential. There we show that the gauge symmetry 
breaking can be achieved consistent with supersymmetry without generating 
pseudo goldstone bosons. Appendix B deals with the threshold corrections 
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to sin 2 6w- 



2. Suppressing Proton Decay in SUSY SO(10) 



The problem with doublet-triplet splitting arises in SU(5) because of the 
tracelessness of its irreducible representations. Suppose the superpotential 
contains the term Xi5' h 24 h 5h and (2A H ) = diag(x, x, x, y, y), then the 3// 
and 3' H higgsinos get a Dirac mass of Airr and the doublets 2 H and 2' H get a 
mass Xiy. We need X±x ~ Mqut and Xiy ~ My/ but this is impossible since, 
by the tracelessness of 24, y = —\x. This can be remedied by introducing 
a singlet superfield, 1#, with the coupling A 2 5' ii l//5// an d = z (or 

equivalently by a bare mass term), but only by tuning the parameters so 
that {-\\ix + \ 2 z) ~ 10- 14 (Aix + \ 2 z). 

In SO(10) such fine adjustment of parameters can be avoided because 
the analogue of the tracelessness condition does not exist. 13 The 24, which is 
the adjoint of SU(5), is contained in the 45 which is the adjoint of SO (10). 
45 is a rank-2 antisymmetric tensor and the VEV of 45// can be brought to 
the canonical form 



which corresponds to the U(5) matrix diag(xi, x 2 , £3, X4, x$). Because this 
is a U(5) rather than an SU(5) matrix its trace need not vanish. One can 
therefore have the VEV of (45//) take the form 



(45//) 



r] <g> diag(xi, x 2 , x 3 , x 4 , x 5 ), 




(3) 
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(45h> = ri® diag(a,a,a,0,0,). (4) 

This is just what is needed to give mass to the SU(3) C - triplet higgs(inos) 
and not the SU(2) L - doublet ones. This is what we call the Dimopoulos- 
Wilczek mechanism. 

There is another group-theoretical explanation for the doublet-triplet 
splitting in 50(10). Under its maximal subgroup SU(2) L xSU(2) R xSU(A) c , 
the standard model singlets of 45 which could acquire GUT-scale VEVs are 
contained in the (1, 1, 15) and (1, 3, 1) multiplets. The 10 of Higgs decom- 
poses as (2, 2, 1) + (1, 1, 6). If only the (1, 1, 15) of 45 acquires a VEV, 
it gives the color triplets of (1, 1, 6) a mass and not the doublets of (2, 2, 
1). If the (1, 3, 1) acquires a VEV, it will supply a super-large mass to the 
doublets - and not to the triplets. (We shall shortly make use of the second 
property to suppress proton decay and at the same time preserve sin 2 6w as 
a prediction.) Such options are not available in SU(5), since SU(5) has no 
intermediate symmetries, even for the sake of classification. 

Consider the following coupling in the superpotential of a SUSY-S'O(IO) 
model: 

W D A10 m 45tfl0 2if . (5) 

One must introduce two 10's of Higgs(ino) fields because with just one the 
term 10#45#10# would vanish by the antisymmetry of the 45. As we noted in 
the introduction, and shall see more clearly later, such a doubling is actually 
a useful thing from the point of view of suppressing higgsino-mediated proton 
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decay. This is another appealing feature of SO (10). 

When the 45^ gets the VEV shown in eq (4) all of the triplet Higgs(ino) 
fields in 10m and IO2H get superlarge masses. The situation can be repre- 
sented schematically as follows 




5m §2H §2H 5ih 

where under 5*0(10) — > SU(5), 10i# = 5i# + 5i#, and 10 2 h = ^>2H + 5 2 h- By 
comparison with the scheme shown in eq. (2) one sees that the 'unprimed 
sector' consists of 5i# and 5 2 #, while the 'primed sector' consists of 5 2 h and 
5i H . One can then identify H = 2 1H and H' = 2i H . 

Now we face the problem of generating mass for the 'extra' doublets which 
reside in 52H + ^>2H — • The simplest possibility is just to introduce into 
the superpotential, W , the term M(10 2 h^0 2 h), with M/M GUT being less 
than - but not much smaller than - one. The resulting threshold correction to 
sm 2 6 w is + 3 " ( ^ z } in ( M GUT /M) « 10" 3 . For proton decay to be suppressed 
it is also necessary that only 10i# but not 10 2 h couple to light quarks and 
leptons with usual strength. All of this (including the absence of a (10i#) 2 
term in W) can be enforced by a global symmetry. For example, one could 
have a symmetry under which 10 2 h — > +IO2H, 10m — > — 10m,45# — > — 45#, 
and 16 j — > HQj (the light families, with j = 1, 2, 3). 

The above appears to us to be the simplest way of achieving weak sup- 
pression of proton decay. (For comparison with SU(5) see the next section.) 
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However, in SO (10), but not in SU(5), it is actually possible to achieve a 
strong suppression in a satisfactory way. To do this we need to give 2 2 // and 
2 2H a superheavy Dirac mass (so as to not mess up sin 2 6 W ) without having 
a superheavy Dirac mass connecting 3 2 // and 3 2 h (which could produce ex- 
cessive proton decay). But this is just a doublet-triplet splitting problem - 
but upside down to the familiar one! Here the doublets but not the triplets 
need a mass term. This will prove to be not doable in SU (5) without tuning, 
but in 5*0(10) it can be done. What is required is another 45 H with a VEV 

(45^) = rj <g> diag(0, 0, 0, a', a') . (7) 

As already noted, this is just as achievable as the VEV given in eq. (4). (See 
section V for the demonstration.) 

There is a slight hitch in that we cannot, because of the antisymmetry of 
45^, simply write down 10^45^10^ to give mass to 2 2 // and 2 2 h- However, 
this can be overcome by introducing an additional 10 of Higgs(ino) fields. 
Consider a superpotential containing 

W D A 1 1H 45 H 1 2H + A' 10 2 ff45' H 10 3 tf 

+ M10 3/ /10 3 // ( 8 ) 

+ v- ; ; /yieae.iom 

with (45//) and (45'#) being given by eqs. (4) and (7). Then the superheavy 
mass matrices of the color-triplet and weak-doublet higgs(inos) are 
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(2i, 2 2 , 2 3 ) 





AW 
-A'a' M 



w 2 ^ 



2 2 
V 2 3y 



(9) 



and 



t 



W 3l \ 



Aa 

(3i, 3 2 , 3 3 ) _ Aa 3 2 ( 10 ) 

m y y 3 3 y 
The doublet matrix is rank- two leaving a single pair of light doublets H = 2i H 
and if' = 2 1H . All triplets get superheavy mass; however, there is no mixing 
between 3i# and 3i# that would permit the diagrams shown in figures 1 or 
3. 

There are several questions to be answered concerning the SO (10) ap- 
proaches to the proton decay problem. (1) Can the VEVs in eqs (4) and (7) 
arise from an actual (super) potential? 14 (2) Can 50(10) be broken all the 
way to SU(3) x SU(2) xU(l) without destabilizing these VEVs? and (3), are 
the threshold corrections in such an S'O(IO) model likely to be small enough 
not to vitiate the successful prediction of sin 2 0^/? We will show in section V 
that the answer to all these questions is 'yes'. But first we will examine the 
possibilities that exist in SU(5) and flipped SU(5). 
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3. Suppressing Proton Decay in SUSY SU(5) 



The only viable method of doublet-triplet splitting in SUSY SU (5) that does 
not involve fine tuning of parameters is the 'missing partner mechanism'. 9 ' 10 
The so-called sliding-singlet mechanism 11 has the problem in SU(5) that 
radiative corrections destroy the gauge hierarchy. 12 For the missing partner 
mechanism in SU (5) one requires (at least) the set 5 + 5 + 50 + 50 + 75 of 
Higgs supermultiplets. In the 50(50) there is a color 3(3) but no weak 2(2). 
Thus the couplings A5#50# (75//) + A'5#50# (75//) give mass to the triplets 
in 5// + 5// but not to the doublets. Schematically, 



/ Q \ 



\2/ 



y other 



/ 



y other 



/3 



\2 



(11) 



5// 50// 50// 5// 

The horizontal solid lines represent superheavy triplet-higgs(ino) masses aris- 
ing from the (75//). As in the cases considered in previous sections, there 
is the question of how to make the 'other' fields in the 50 + 50 superheavy. 
(They contribute to the RGE at one loop the same as a pair of weak dou- 
blets.) If one wanted a 'strong suppression' of proton decay, it would require 
giving superlarge masses to all the 'other' fields in 50 + 50 but having the 
3(50//) and 3(50//) not be connected by a large Dirac mass term. There is 
no analogue of the missing partner mechanism that would accomplish this 
in a natural way. It could only be done by fine-tuning. For example, two 
different representations (1// and 24//, or 24// and 75//) could couple 50 to 
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50 and be relatively adjusted to give light mass only to the 3(50#) + 3(50#). 
However, we have foresworn fine-tuning. 

It is possible to achieve a natural weak suppression of higgsino-mediated 
proton decay in SU(5) by introducing an explicit mass term M(50#50#) into 
the superpotential and having M/Mqut be of order but somewhat smaller 
than unity. This works as well as the weak suppression mechanism in SO (10) 
discussed in section II. However, in SU(5) there is the necessity of introducing 
the somewhat exotic high rank representations 50, 50 and 75, whereas in 
SO (10) only the low rank representations 10, 45, and 54 are required. If one 
were willing to live with multiple fine-tunings of parameters one could do 
with just doubling the higgs sector in SU (5) to 5 + 5 + 5' + 5' as discussed 
briefly in the introduction. With two fine-tunings one could make all the 
triplets heavy, and achieve weak suppression of proton decay. With a third 
fine-tuning one could give mass to the extra pair of doublets and yet achieve 
strong suppression of proton decay. There are papers in the literature that 
take this approach. 7 ' 8 

4. Suppressing Proton Decay in Flipped SU(5) 

As is well-known, the missing partner mechanism works much more eco- 
nomically in flipped SU(5) 15 than in ordinary SU(5). 16 One requires for the 
mechanism the SU(5) x U(l) representations 5^ 2 + 5 2 H + 10# + 10^/, and 
the superpotential couplings A5# lO^lOjj + A^IO^TO^ 1 . The 10 h (10h) 
contains a color 3(3) but no color-singlet, weak-doublet components. The 
10h and 10 h get VEVs that break SU(5) x U(l) down to Gs and also mate 
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the triplet higgs(inos) in the 5h + 5# with these in the IOh + 10 h leaving 
the doublets in 5# + 5h light. Schematically, 



(') 


[ : ) 


[ ' ) 


[■) 


V 2 y 


y other J 


\ other J 


v 2 y 


5 H 




VOh 


5h 



(12) 



Another beautiful feature of flipped SU(5) is that there is no necessity to 
do anything else to give mass to the 'other' fields in the IOh + IOh- they 
are all disposed of by the (super) Higgs mechanism! They are either eaten 
or become superheavy with the gauge/gaugino particles. Thus, in flipped 
SU(5) one can strongly suppress higgsino-mediated proton decay without 
any fine-tuning and without leaving any 'extra' split multiplets lighter than 
M G ut- We found this to be impossible in ordinary SU(5). However, there 
is one major drawback: the group of flipped SU(5) is really SU(5) x U(l) 
and so real unification of gauge couplings is not achieved. One has therefore 
lost, or rather never had, the unification of gauge couplings as a prediction. 

5. A More Detailed Examination of SO(10) 

In section II certain ideas were discussed for solving the doublet-triplet split- 
ting problem and for suppressing higgsino-mediated proton decay in S'O(IO) 
that made essential use of specific patterns of VEVs, in particular those 
shown in eqs. (4) and (7). The question arises whether such VEVs are nat- 
ural. In ref. 14 Srednicki wrote down a superpotential for a 45 and a 54 of 
Higgs that has both of these forms as possible supersymmetric minima. Let 
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us denote the 45 and 54 by A and S respectively. Then the most general 
S'O(10)-invariant superpotential involving just these fields has the form 

W(A, S) = m x A 2 + m 2 S 2 + AiS 3 + X 2 A 2 S . (13) 
The equations for a supersymmetric minimum are 



= F A = 2{ mi + X 2 S)A 

= F S = (2m 2 S + \ 2 A 2 + Z^S 2 ) . (14) 



1 



Suppose we choose (S) — I .". " 1 <S> diag(s, s, s, — |s, — |s) (S is a traceless 



l I " ") ") 2"' 2' 

-1 



rank-two symmetric tensor of 5*0(10)) and (A) = ( ^ ^ J ®diag(a, a, a, 6, 6). 



Then equation (14) gives two equations: 

(mi-§A 2 s)6 = 0, 

(15) 

(mi + A2«)a = 0. 
Either a or 6 or both must therefore vanish (if s ^ 0). There are therefore 
three possible solutions. (1) b = 0, s = —mi/\ 2 ; (2) a = 0, s = ^mi/\ 2 ; (3) 
a = b = 0. 

For doublet-triplet splitting and weak suppression of proton decay we 
need only the solution (1). For strong suppression we need (at least) two ad- 
joints, one with VEV corresponding to solution (1) and the other to solution 
(2). (See eqs. (4) and (7)). We will examine this latter more complicated 
case in greater detail. If the required pattern of VEVs can be achieved in a 
realistic model for that case, then a fortiori the simpler requirements for weak 
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suppression can be achieved also. The main issues are whether the VEVs 
in eqs. (4) and (7) can be achieved, the group SO (10) broken completely 
to SU(3) x SU(2) x U(l), and goldstone particles avoided. (The issue of 
threshold corrections to sin 2 9w is dealt with in Appendix B.) 

To begin with we double the superpotential shown in equation (13). That 
is, we have two 45's, denoted A and A', and two 54's, denoted S and S', with 
superpotential 

W(A,S;A',S') = m 1 A 2 + m 2 S 2 + X l S 3 + X 2 A 2 S 

(16) 

+ m! x A 12 + m' 2 S' 2 + \[S' 3 + \' 2 A' 2 S'. 
We aim to have the (A, S) sector have VEVs in solution (1), and the (A', S') 
sector have VEVs in solution (2). 

The superpotential in (16) is certainly not enough because, for one thing, 
nothing determines the relative alignment of the VEVs of the two sectors, 
and so there are goldstone modes corresponding to a continuous degeneracy 
whereby the sectors are rotated in 5O(10)-space with respect to each other. 

A second problem with eq. (16) is that (A) and (S) break 5*0(10) — > 
[SU(3) x 1/(1)] x 50(4), while (A') and (S') break 50(10) -> 50(6) x 
[SU(2) x U(l) x U(l)]. Altogether, then, the unbroken group is SU(3) x 
SU(2) x 17(1) x 17(1) = rank 5. 

To break all the way to the standard model further Higgs fields are needed. 
(They are needed for right-handed neutrino masses in any case.) The simplest 
choices are 16 + 16 or 126 + 126. Let us call these C + C where C = 16 or 
126. One can write a superpotential that gives C + C VEVs which break 
517(10) -> SU(5). Together with A, 5, A', and S' this will complete the 
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breaking of SO (10) to Gs and allow ur masses. 

At this point a further somewhat subtle technical problem arises. There 
are certain generators of SO (10) that are broken both by C + C and by 
the adjoints A and A', but not by the symmetric tensors S + S', specifically 
the generators in SO(Q)/SU(3) x U(l) and SO(A)/SU(2) x U(l) (where 
SO (10) D 5*0(6) x SO (A)). Thus to avoid residual goldstone bosons there 
must be coupling between the C + C sector and the adjoints A and A'. 
The technical problem is that the direct couplings CAC and CA'C would 
destabilize the desired VEVs of A and A'. In particular, all the 'diagonal' 
components of (A) and (A') (written as U (5) matrices, that is) become non- 
vanishing. [This is because the VEVs of C + C that break SO (10) — > SU(5) 
couple in CAC to the S , f/(5)-singlet combination (3a + 26), where (A) = 



(3a + 2b) 2 in the ordinary scalar potential, which in turn leads to cross-terms 
of the form ab, destabilizing the solution a ^ and 6 = 0. The same thing 
would happen to A'.] 

There are various solutions to this technical difficulty The one we shall 
study here involves the introduction of a third adjoint, which we will denote 
A", that serves as an intermediary between the C + C sector and A and A'. 

The part of the superpotential that does the complete breaking to the 
standard model and obviates all difficulties is given in full by 




x diag(a, a, a, b, b). This leads to a mass term proportional to 
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W = mi A 2 + m 2 S 2 + \ 1 S 3 + \ 2 A 2 S 
+ m^A' 2 + m' 2 S' 2 + A'^' 3 + \' 2 A' 2 S' 

(17) 

+ m'(A" 2 + m' 2 'CC + \%CA"C 
+ XAA'A". 

There are three sectors, (A, S), (A' , S'), and (A",C + C), that are coupled 
together only by the last term XAA'A". 

The term X'^CA'C does serve to give A" a VEV in the S'C/(5)-singlet 
direction: 

(A") = ( J ) ® diag(a", a", a", a", a"). (18) 

But this does not destabilize the VEVs of A and A' which are assumed to 
be of the forms given in eqs. (4) and (7). This is easily seen by examining 
the XAA'A" term, which is the only thing linking A" to A and A'. Consider 
the Fa = equation. (Fa)^ is an antisymmetric tensor to which XAA'A" 
contributes X(A'A"Y ab \ which vanishes when the values of (A') and (A") 
given in eqs (7) and (18) are substituted. Similarly, at the desired minimum, 
XAA'A" gives no contribution to the Fa 1 and Fa" equations. In other words 
it can be neglected in doing the minimization! However, it does contribute 
to the Higgs(ino) masses, and, indeed, removes all of the possible goldstone 
modes discussed above. 

In Appendix A we present details of the minimization of the superpo- 
tential, eq. (17), assuming C = 16. There we show explicitly that SO (10) 
may be completely broken to the standard model (without breaking SUSY), 
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uneaten goldstone bosons avoided, and VEVs of the desired form achieved. 
The masses of the various Higgs (super) multiplets enumerated in Appendix 
A will be used to estimate threshold corrections in the model. 

It is conceivable that the superpotential of eq. (17) is the most general one 
compatible with some discrete symmetries, although due to the supersym- 
metric non-renormalization theorem one is not obliged to write all possible 
terms. 

We have found other superpotentials and sets of Higgs fields that allow 
us to achieve the desired VEVs in a consistent and realistic way. We have 
presented eq. (17) as being algebraically simple to analyze. It should also be 
noted that the implementation of weak suppression of HMPD, where only a 
single 45 h is needed with VEV of the form in eq. (4), is a simpler task and 
fewer fields are required. We have not tried to find the absolutely minimal 
scheme. 

At this point we wish to make an aside. If one is willing to give up 
sin 9\y as a prediction, there is a much simpler way to simultaneously achieve 
doublet-triplet splitting without fine-tuning and a strong suppression of hig- 
gsino mediated proton decay. All we need is then one 45 H of Higgs superfield 
with its VEV as given in eq. (4). Suppose the relevant superpotential term 
for doublet-triplet splitting is just A10i#45#102# as in eq. (5). This will 
make the color triplets heavy, but one is left with 2 pairs of light doublets. 
Now, if the 45# does not couple to the sector that breaks 5*0(10) — > SU(5) 
(via C + C superpotential), then in addition to the extra pair of doublets, one 
will have a {(3, 1, |) + H.C.} goldstone super-multiplet which remains light. 
(These are the goldstones corresponding to SO (6) /SU(3) x £7(1) mentioned 
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earlier.) The combined effect of having an extra pair of light higgs doublets 
and the charge | higgs (super) fields is to alter sin 2 9 W prediction to ~ 0.215 
at one loop. The unification scale also comes down by an order of magnitude 
or so. If one 'fixes' these features by relying on particle thresholds, such a 
scenario may not be inconsistent. This scenario can be tested by directly 
searching for the (3, 1, |) - higgs and higgsino particles. This situation is 
somewhat analogous to the case studied in ref. 17. We do not advocate 
this scenario here, since our aim is to preserve the successful unification of 
couplings as a prediction. 

Returning to the superpotential in eq. (17), it might be imagined that 
with 3(45) + 2(54) + 16 + 16 the threshold corrections might be fairly large; 
large enough, perhaps to vitiate the successful 'prediction' of sin 2 # w . Actu- 
ally, this is not the case, especially if one assumes SO(10) breaks in two stages 
to the standard model: SO(10) -> SU(5) -> SU(3) x SU{2) x U(l). SO(10) 
is broken to SU (5) by the VEVs of C, C, and A" at a scale M 10 , while SU (5) 
is broken down to SU(3) x SU{2) x 17(1) by the VEVs of A, S, A', and S' at 
a scale M 5 . The masses of particles will be of the form aM 10 + @M 5 . In the 
limit that (3M 5 /aM 10 — > for a given multiplet its one-loop threshold cor- 
rections to sin 2 9 W will vanish since it will become a complete and degenerate 
SU(5) multiplet. Thus threshold corrections of complete S'C/(5)-multiplets 
go as £n(aM 10 + (3M b )/{aM w ) = (3M 5 /aM w for M 5 < M 10 . Thus if M 10 
is assumed to be somewhat larger than M 5 , the GUT-scale threshold correc- 
tions to sin 2 6 W are substantially reduced. These will be discussed explicitly 
in Appendix B, where it is found that the uncertainties in sin 2 9w due to 
superheavy thresholds is typically in the range of 3 x 1CT 3 to 1CT 2 . 
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6. Conclusions 



If one seeks a supersymmetric grand unified model in which proton decay 
mediated by color-triplet higgsino is strongly suppressed through a mecha- 
nism based on symmetry, in which there is no fine-tuning of parameters, and 
in which the remarkable prediction of sin 2 9w is maintained as a prediction, 
then it seems that one must turn to SO (10). On the other hand, a weak sup- 
pression due not to symmetry but to the smallness of a parameter is achiev- 
able in both SU(5) and SO (10) without either fine-tuning or sacrificing the 
sin 2 9 W prediction, though we believe SO (10) allows the more economical so- 
lution. The SU (5) solution, being based on the 'missing-partner mechanism', 
requires the introduction of Higgs(ino) multiplets in 50 + 50 + 75 (which are 
four and five index tensors), whereas the 50(10) solution requires only the 
usual (rank one and two) tensors 45, and 54, and the spinors 16 + 16. 

The advantage of SO (10) is due to the possibility of exploiting the elegant 
Dimopoulos-Wilczek mechanism of doublet-triplet splitting. We have studied 
that mechanism and found that it can be implemented in a fully realistic way. 

In our view these results constitute yet another argument in favor of 
50(10). It is already well known that 50(10) has the advantage over SU(5) 
of allowing .R-parity to be a gauge symmetry (that is because Higgs fields 
are in tensor representations and matter fields are in spinor representations). 
And, of course, 5O(10) achieves greater unification of quarks and leptons 
and requires the existence of right-handed neutrinos. 

In any event, we have shown that higgsino- mediated proton decay is not 
a serious difficulty of supersymmetric grand unification as there are quite 
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simple and natural means to suppress it without undercutting the main suc- 
cess of those models. If the suppression is of the 'weak' type then there are 
grounds to hope to see proton decay in super Kamiokande. 
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Appendix A 



In this Appendix, we give details of the minimization of the superpotential 
of eq. (17). We shall see explicitly that (a) SO (10) breaks completely to the 
standard model in the supersymmetric limit, (b) the desired forms of the 
VEVs of A and A' are achieved, and (c) there are no unwanted pseudo gold- 
stone modes which could potentially ruin the successful sin 2 9w prediction. 

A 2 in eq. (17) denotes Tr(A 2 ), AAA" denotes Tr(AA'A") etc. We shall 
confine ourselves to the case where C + C = 16 + 16. The term CA"C is 
explicitly written down as Ca^A^C '/4, where a a p are the generators of 
SO(10) algebra. 18 

The VEVs for the fields are chosen as 

(A) = 7] <S> diag(a, a, a, 0, 0) 
(A') = 7] <g> diag(0, 0, 0, a', a') 
(A") = V ® (a", a", a", a", a") 

(A.l) 

(S) = / ® diag(s, s, s, -§s, -§s) 

(S') = /(8>diag(s',s / ,s / ,-|s',-|s / ) 
(C) = (C) = c 

where rj = ^ ^ J ^ and I — ^ J ^ ^ . The equality (C) = follows 
from the requirement of D-flatness. The vanishing of the F-terms lead to the 
following conditions, corresponding to A, S, A', S', A", C fields respectively: 
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mi + X 2 s 

\ 2 1 \ 2 



m', — -Aos' 

1 2 2 

3 1 

/ / \ / /2 i \ / /2 

m 9 s — -AiS H — A 9 a 

2 4 i 5 2 

\" 

^//^// | A 2 2 



= m 



2 2 -2 



An a . 



(A.2) 



Since SUSY is unbroken, it is sufficient to investigate the Higgsino mass 
spectrum. The multiplets which transform as {(3,1,2/3) + H.C.} under 
SU(3)c x SU(2) L x U(l)y have the following mass matrix: 



Mi 







2Aa" 











2Aa" -10A' 2 6' 



(A.3) 



2Aa 

-2Aa -Am'[ -A' 2 'c 

V -A 2 'c 2A 2 'a" / 

This matrix has one zero eigenvalue by virtue of eq. (A.2). All the other 3 
states become massive. 

The mass matrix corresponding to {(1, 1, +1) + H.C.} is 
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/ lOXob 2Xa" -2Xa' \ 



2\a" 















(A.4) 



-2Xa' -Am'l -A' 2 'c 

V o o -a 2 ' c 2\y j 

Again, M. 2 has one zero eigenvalue (using (A. 2)) and three non-zero eigneval- 
ues. 

The mass matrix for {(3, 2, — |) + H.C.} is given by 



( 5A 2 s V2\ 2 a 
\f2\ 2 a 2m 2 - |Ais 



-5A' 9 s' 






V2X' 2 a' 



-Xa' ^ 


— Xa 



(A.5) 



V2X' 2 a' 2m 2 -|A , 1 s' 

V -Xa' -Xa -4m'/ / 

Using (A. 2) one sees that has one of its eigenvalues equal to zero, while 
the rest are all nonzero. 

The corresponding matrix for {(3, 2, |) + H.C.} is given by 
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Ma = 



5X 2 s 



\[2X 2 a 2m 2 — f Ai« 



-2i\a" 



iXa' 













2i\a" 



-5X' 2 b' V2X' 2 a' 

V2X' 2 a' 2m' 2 - fA'^' 

— iXa 




—iXa' 


iXa 


-Am'l 
X'ic 










A' 2 'c 



2X'ia" j 
(A.6) 



This has one zero and five nonzero eigenvalues. 

The goldstone modes in .Mi, -M 2 , .M3, -M 4 when combined with the zero 
mass mode corresponding to the phase of (C + C')-singlet add up to the 
33 mass-less modes needed for the symmetry breaking S'O(IO) — > SU(3) x 
SU(2) x U(l). All the remaining fields become massive. Their spectrum 
looks as follows. From the (A, S) sector, we have 



{(6,1,--) + H.C.) = 2m 2 + 6X lS 
{(1,3, ±1); (1,3,0)} = 2m 2 -9X lS 



{(1,3,0) + (1,1,0)} = 10A 2 s 

I 2V2X 2 a \ 



{(8,1,0)} = 



\ 2v / 2A 2 a 2m 2 + 6X1S J 
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{(1,1,0)} = 
From the (A', S f ) sector, one finds 



\ 



y — 4?A 2 a 2m 2 — 3AiS j 



(A.7) 



{(6,l,--) + tf.C.; (8,1,0)} = 
{(8,1,0); (1,1,0)} = 



2m 2 + 6Xis' 



-4A' 2 s' 



{(1,3, ±1)} = 2m , 2 -9A' 1 s / 



{(1,3,0)} 



{(1,1,0)} = 



-2V2\' 2 a' 



-2^2\' 2 a' 2m' 2 -9\' 1 s' 
( n /MV„' \ 





f2A\l „l 



. (A.8) 



^ A' 2 a' 2m' 2 - 3A' lS ' ) 
Finally, from the (A", C + C) sector, one finds 



{(8, 1,0) + (1,3,0) + (1,1,0)} 



-Am" 



{(S,1,~) + H.C; (1,2,^) + H.C.} = 2m' 2 ' + 3A' 2 'a" 



{(1,1,0)} 



-4m'/ 



\ 



-J|A 2 'c m' 2 ' 



|A 2 a" 



(A.9) 



Now to reduce threshold corrections somewhat (and simplify calculations) 
we assume the scale of 5*0(10) breaking, M 10 , is somewhat greater than 
the scale of SU(5) breaking, M 5 . (This means we assume m'(,m 2 ,c,a" ^> 
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mi,m 2 ,m' l ,m' 2 ,a,a',s,s'.) As explained in the text, the multiplets which 
have mass 0(M W ) will give contributions to the threshold corrections sup- 
pressed by M 5 /Mi . Thus in Appendix B we will only need the masses of 
particles which are 0(M 5 ). For sets of particles with the same Gs quantum 
numbers we only will need to know the products of their masses. These are 
listed below. 



{(3, 2, |) + h.c.} = {±)\ 2 \ 2 a 2 a' 2 /ss' (2 states) 
{(3,2,-§) + /i.c.} = (±)(\ 2 \ 2 \' 2 /\"(a 2 a'a"/ss'c 2 ) 

\2aa!(a + a') + 25(a's 2 + as' 2 )} (3 states) 

{(6,l,-§) + /i.c.} = (f)X 1 s+(l)X 2 a 2 /s (1 state) 

{(1, 3, ±1); (1,3,0)} = -(f)A lS + (|)A 2 a 2 /s (1 state) 

{(1,3,0)} = 10A 2 s (1 state) 

{(8,1,0)} = 8\ 2 2 a 2 (2 states) 

{(6,1,-1) + h.c; (8,1,0)} = (f )A / 1 s / - (§)A 2 a /2 /s' (1 state) 

{(8,1,0)} = -A\' 2 s' (1 state) 

{(1,3,±1)} = -(f )Ais' - (f )A 2 a /2 /s' (1 state) 

{(1,3,0)} = 8A' 2 V 2 (2 states) 

(A. 10) 



We have expressed these in terms of the VEVs and Yukawa couplings and 
eliminated the mass parameters, m\,m 2 , etc., using eqs. (A. 2). 
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Appendix B 



Here we use the results of Appendix A to compute the threshold corrections 
to sin 2 6w coming from superheavy fields. We look at the full model with 
three 45's (A, A', A") and two 54's (S, S'), and we assume 50(10) breaks to 
SU(5) at a scale (M w ) which is higher than the scale at which SU (5) breaks 
to the standard model (M 5 ). 

It will prove convenient to define the parameters M = A2S, x = a/s,x' = 
a'/s', y = A1/A2, y' = A^/A^, and z = X' 2 s' '/X2S. The correction to sin 2 9w is 
given at one loop by 



A sin 2 6 W {M Z ) = ^± [5b{ - \2V 2 + 7b> 3 ] In M j (B.l) 

where the sum is taken over multiplets of SU(3) x SU(2) x U(l). From eqs. 
(A. 10) one obtains 



Asin 2 MM z ) = ^{-211n(^xV 2 zM 2 ) 

+ 31n(tM 3 ) +51m((f y+ \x 2 )M) 



(B.2) 



+ (-30 - 24) ln((-f y + |x 2 )M) - 24 ln(10M) 

+ 21 ln(8a; 2 M 2 ) + (51 + 21) ln((f y' - \x l2 )zM) 

+ 21 \n(4zM) - 30 ln((-f y' - \x l2 )zM) 
- 241n(8x ,2 z 2 M 2 )}. 

(Here tM 3 is defined to be equal to the complicated expression on the right- 
hand side of the second equation of (A. 10). t is of order M 5 /M w . However, 



30 



the coefficient of this term is mercifully small, so that its effect is negligible. 
We have not displayed the threshold effects due to the doublet and triplet 
fileds of eqs. (9)-(10) as they are negligible.) All logarithms are understood to 
have an absolute value in their arguments. We have used 63(3) — \, 63(6) = 
|, 6s(8) = 3, b 2 {2) = I, 6 2 (3) = 2, ^(f) = |(y/2) 2 . 
Collecting terms, 

Asm 2 6 W (M Z ) = ^ 1 {(181n5-331n2) 
+ 31nt-61n^-451n(x' 2 ) 

+ 511n(fy + §:r 2 )-541n(fy-§a; 2 ) ■ ( R3 ) 

+ 72 ln(f y> - \x 12 ) - 30 ln(f y> + fx' 2 ) 
- 31nM} 

The term — 31nM is due to the (3,2, — |) + h.c. that get eaten when SU(5) 
breaks and is present also in the minimal SU (5) model. If it were not for this 
eating, all the superheavy multiplets would be complete SU (5) multiplets and 
the dependence on M would disappear. (It is only splitting within multiplets 
that contribute at one loop to sin 2 Q w . Since all masses are scaled by M, M 
drops out in the ratios.) In the first term, 18 In 5 — 331n2 = 6.1. All of the 
logarithms can have either sign. There are five potentially large terms with 
coefficients averaging about 50. If we assume the logarithms are of order one 
with arbitrary signs then the typical threshold correction to sin 2 9w would be 
expected to be about ±(a/307r)(10 2 ) ~ ±10 -2 . This is to be contrasted with 
the effect of a pair of extra light Higgs doublets of A sin 2 9 W ~ 2.5 x 10 -2 . It 
should also be compared to the theoretical uncertainties in sin 2 9w in minimal 
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SUSY SU(5), referred to in the opening paragraph of this paper, of about 
0.36 x 1(T 2 . 

The expression in (B.2) simplifies considerably if we assume y < j^x 2 , y' < 
±x 12 . Then 



Since 6 In 2 — 21 In 5 = —29.4 one expects the threshold correction to be 
negative and about —0.3 x 10 -2 in this limit. We mention this limit since 
it is a special solution of the superpotential of eq. (17), corresponding to 
setting the parameters Ai, \[ to zero. We note that this limit can be reached 
naturally without giving rise to any pseudo goldstones (see Appendix A). 

We conclude that the threshold correction uncertainties to sin 2 9 W in the 
kind of 50(10) model we are discussing are likely to be a few times larger 
than the total theoretical uncertainty in sin 2 6V in minimal SUSY SU(5), 
but several times smaller than the effect on sin 2 9 W of an extra pair of light 
Higgs doublets. We should emphasize that if one is satisfied only to suppress 
Higgsino mediated proton decay weakly, a much smaller Higgs sector may 
be adequate, with correspondingly smaller threshold corrections. 
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(B.4) 



6 In z + (61n2-21hi5)} 
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Figure Captions 

1. A diagram that gives a dimension 5 baryon- number- violating opera- 
tor. The arrows indicate the direction a left-handed chiral superfield is 
flowing. Chirality shows this to be an F term and hence suppressed by 
(Mass)^ 1 . Chirality also implies that there must be a mass insertion 
(denoted by M) on the colored higgsino line. The suppression is thus 
M/Mq UT which is naturally of order M^j T in most models. 

2. A diagram without the chirality-flipping mass insertion of Fig. 1 and 
thus representing a D term. Such a term is effectively of dimension 
6 and suppressed by Mq^ t . It corresponds to colored higgs boson 
exchange. 

3. In the 'weak suppression' scheme discussed in section I (see eq. (2)) the 
primed and unprimed Higgs(ino) sectors are connected by a Dirac mass, 
M. This coupling of the two sectors allows a dimension-5, baryon- 
number-violating operator to arise as shown here. 
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